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I discuss the pt distributions for high-multiplicity events 
originating from semi-classical variation of the gluon density 
of the proton. The multiplicity distribution measures the cur- 
vature of the effective action for the small- a; gluon fields. For 
pA collisions at the RHIC and LHC colliders, semi-classically 
the multiplicity distribution reflects the distribution of sat- 
uration momenta of the proton but not that of the nucleus. 
The average transverse momentum in the central region grows 
with dN/dy, while the pt distribution of leading hadrons in 
the proton fragmentation region should depend less on the 
multiplicity in the central region. 

High-energy hadronic scattering requires understand- 
ing of the non-abelian gauge fields of hadrons at small x, 
which is the fractional light-cone energy carried by the 
quanta of the fields. At very high energy, logl/x ^> 1, 
the gluon fields in a hadron become very strong, cor- 
responding to high gluon density. This is where one 
expects that cross sections become comparable to the 
geometric size of the hadron and where the unitarity 
limit is reached. A perturbative QCD based mecha- 
nism for unitarization of cross sections is provided by 
gluon saturation effects Q . A semi-classical approach to 
gluon saturation and QCD at high energy was developed 
in [^|J|^,||,^,^U|] and will be applied here to high energy 
proton- nucleus collisions. 

Large nuclei facilitate the study of gluon saturation ef- 
fects because the gluon density per unit transverse area 
is larger than in a proton. The scale associated with the 
high gluon density, the saturation scale Q s , grows with 
energy and atomic number A. At a resolution less than 
Q 2 , the color field carries large occupation numbers, of 
order of the inverse QCD coupling constant, l/a s . Thus, 
the nuclear wave function at Q 2 < Q 2 resembles a Bose 
"condensate" . The local color charge density in the trans- 
verse plane is a stochastic variable which eventually has 
to be averaged over, see below. Also, by Lorentz time di- 
lation the large x gluons evolve slowly and so for the 
small- x gluons they appear as a "frozen" source near 
the light cone. Therefore, the high gluon density state 
of QCD at Q 2 < Q 2 is called a "Color Glass Conden- 
sate" §. 

The small- x gluon fields of hadrons or nuclei can be 
treated semi-classically because the occupation numbers 
parametrically are of order l/a s S> 1. One can inte- 
grate out the fields at large x whose dynamics is "frozen" , 
thereby generating an effective action for the small- a; glu- 
ons iig|i0§, 

S[p] = [ dyd 2 x t V[p] , (1) 



with y being the rapidity. The large-x gluons effectively 
act as a source of color charge in the Yang-Mills equations 
of motion for the small- x fields; p denotes the color charge 
density per unit transverse area d 2 Xt, and rapidity dy. 
This is a stochastic variable which eventually has to be 
integrated over, 

{0) = J Vpex P (-S[p])0[p] , (2) 

where 0[p] denotes some observable which is a functional 
of p. In their original papers McLerran and Venu- 
gopalan suggested the simplest possible effective poten- 
tial 

V 2 [p;p = —2- ■ 3 
M 

This is the lowest-dimensional operator. In (H|), p 2 is a 
real number, and physically is simply the mean square 
fluctuation of the density of "hard" , large- a; gluons in 
the source. It determines the curvature of the effective 
potential at the minimum. 





V[p] 






1 

i 


i \ / 

^2 \ / 

\ / 

i / 



R ft 



FIG. 1. Double- well potential for the color charge density 
per unit area. 

In principle, one might try other choices for the effec- 
tive potential V[p], c.f. [p|JlO|, as for example the double 
well potential depicted in Fig. [j], which looks just like 
the effective potential for a first-order phase transition in 
thermal SU (3) gauge theory near T c [jllj . The curvature 
of the potential in the first well at p\ is 1/pi, and that 
in the second well is l//i2 ■ Semiclassically, one considers 
only small fluctuations about the states which minimize 
the action. Thus, for this potential the averaging over p 
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from eq. (0) turns int 



J 



Wj / Vp exp ( - / dy d 2 x t V 2 [p - \ 0[p] . 

z=l,2 \ J J 

(4) 

The probabilities for the configurations pi or p 2 are de- 
termined by the classical action: 



exp(-5[pi]) 



exp(-5[pi]) + exp(-S[p 2 ]) 



(5) 



Note that if the effective potential is quadratic for small 
fluctuations about each of its minima, renormalization- 
group evolution presumably preserves that functional 
form HHim . That is, imposing the potential from Fig. |l| 
as the initial condition (at small rapidity) , it will keep its 
shape at larger rapidity, and the ratio of the curvatures 
at the two minima stays the same. 

One can evidently generalize the above to effective po- 
tentials with several classical minima. One could also 
imagine that there is a continuum of nearly degenerate 
states near a classical minimum. For example, at a tri- 
critical point the lowest-dimensional operator appearing 
in the action is ~ trp 6 . The state of lowest action is that 
for p = but the potential about that classical state is 
flat, and so states in its vicinity are nearly degenerate. 
This type of potential might not be preserved by the RG 
evolution in rapidity, though. 

The distribution w(p 2 ) of correlation lengths deter- 
mines the multiplicity distribution in hadronic collisions. 
Here, I do not attempt to compute that distribution. 
Rather, given a nonvanishing probability for a class of 
events with higher multiplicity dN/dy than the average, 
I discuss the origin of such events and how the pt dis- 
tribution in that class of events is different from the av- 
erage. Nevertheless, it will be very interesting of course 
to analyze the multiplicity distribution experimentally; 
for example, at a tricritical point events with larger than 
average multiplicity could be due to semi-classical fluc- 
tuations of the color charge density. A potential with 
several discrete local minima as in Fig. ^ would lead to 
multiple peaks in the multiplicity distribution, unless the 
curvature radii of the potential are the same. For finite 
temperatures, the effect of a second local minimum in the 
effective potential on multiplicity and momentum distri- 
butions of produced particles has been discussed previ- 
ously using the hydrodynamical model |l2|| . 

Henceforth, I consider "pA" collisions, i.e. collisions of 
a dilute projectile with a dense target. For simplicity, p 2 



1 The operator O is assumed to be defined such that tadpole 
contributions are subtracted; for example, for the two-point 



for the nucleus is chosen to be constant over the trans- 
verse plane. For large nuclei (A ~ 200) this should not 
be a bad approximation. Realistic transverse density dis- 
tributions can be treated numerically Moreover, ex- 
perimental triggers in the fragmentation region of the nu- 
cleus, for example a trigger on the number of knocked-out 
neutrons, could improve the selection of central events. 
This might be particularly important for composite pro- 
jectiles (small nuclei) to ensure that all nucleons have hit 
the large nuclear target. 

For u pA v collisions an analytical solution for the ra- 
diation field in the forward light-cone can be found to 
first order in the coupling to the proton field, but to all 
orders in the density of the nucleus (see below). In prin- 
ciple, there is a contribution to large multiplicity events 
from coupling the radiation field to second (and higher) 
order to the proton source, down by additional powers 
of a s . That corresponds to attaching additional non- 
collinear hard gluon lines to the proton source (collincar 
emissions are already resummed in eq. ( |l5|) through the 
DGLAP logarithm), see e.g. H]. As already mentioned 
above, for the simplest quadratic effective potential this is 
in fact the only contribution to large-multiplicity events 
since the effective potential exhibits only one classical 
minimum (i.e. the distribution of saturation momenta 
is a 5- function) . In general, if the potential starts out 
with some higher-dimensional operator, the distribution 
of saturation momenta is less sharply peaked, or it ex- 
hibits several peaks as for the double well potential above, 
and so high-multiplicity events can occur semi-classically. 
Their contribution dominates when exp(— AS'[p]) > a s , 
where A5[p] is the action relative to that at the global 
minimum. This condition determines the "cut-off" of the 
semi-classical multiplicity distribution. 

In light-cone gauge, the classical fields of the sources 
before the collision, which occurs at the tip of the light- 
cone, are given by the single-hadron (nucleus) solutions 

Af <2 = , 

4,2 = -Utfl^&UliiXA.) . (6) 
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The Lq^'s are rotation matrices in color space, 

Ui, 2 (x t ) = exp (-ig$x,2(x t )) . (7) 

The gauge potentials satisfy the two-dimensional 
Poisson equation 



-V 2 $i, 2 (x t ) = gpi, 2 (x t ) 



(8) 



function (p(x t ,y)p(z t ,y')) 
9^ p + Pi 



Pi . One can then shift variables 



The average over the color charge density configurations 
pi and P2 with a Gaussian weight as in eq. (|^) thus cor- 
responds to an average over the rotations in color space, 
C/1,2. 

The fields (^) serve as the initial conditions for solving 
the Yang-Mills equations of motion in the forward light 
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cone; they were solved in |L5| for the case of asymmetric 
collisions where the classical field of one of the colliding 
sources is much stronger than that of the other source. 
Chosing the gauge x + A~ + x~A + = 0, the radiation field 
in the forward light cone region is given by 

A\T,x t ) = U{x t ) (p l (T,x t ) + -A U\x t ) , 
A ± {r,x t ) = ±x ± U{x t )!3{T,x t )U\x t ) . (9) 



at central rapidity, for gluons produced with transverse 



momenta Q p (y) 



< 



Pt 



< 



Qa(v) the proton projectile is 



dilute while the nuclear target is dense. Performing the 
integrals over color charge density configurations with 
a quadratic potential for fixed saturation momenta Qa 
and Q p , the transverse momentum distributions in the 
two regimes are given by |~5]]a 



dN l Q 2 A (y)Q 2 (y). P 2 



cPbcPpt dy g 2 pf 



Q A (y) 



Here, t = \/2x + x~ denotes proper time and x± is the 
transverse coordinate. The £/'s are rotation matrices in 
color space, to be specified shortly. At asymptotic times, 
r — * oo, the fields (3 and (3 l are given by superpositions 
of plane wave solutions, 



(Pt > Qa) 



dN 



i Qliv) , P 2 t 



d 2 bd 2 p t dy g 2 pf 



log- 



(14) 



(15) 



(Q P <Pt< Qa) ■ 



0{t 



oo, x t ) = 
d 2 Pt 



1 



(2tt) 2 72^3 
oo,x t ) = 
d 2 p t 1 



{ai(Pt)( 



-iPfXt—lCJT 



+ C.C.} , 



^ l p\ 



{tt2(Pt)( 



ipt -Xt—IUJT 



Again, some proportionality constants which are inessen- 
tial for the present discussion have been dropped 
(10) ( see 101 )■ The logarithm in eq. ( [l5| ) arises from resum- 
ming collinear gluon emissions from the proton source 
at the scale p\ jl(|, i.e. from DGLAP evolution. For 
Qa 3> Qp — Aqcd the dominant contribution to the mul- 
+ c.c.} . (11) tiplicity density per unit rapidity is from the region (ITF 



The number distribution of produced gluons at rapidity 
y and transverse momentum p t is given by 



dN 



f P (y) log 



dN 



d 2 Pt dy (2tt) 2 



tr {\ai(pt)\* + \02(pt)Y 



(12) 



d 2 bdy g 2 ^ pyy ' b Q 2 p (y) ' 
The contribution from the high-p t region (|l4j) is 
dN 1 



For pA collisions, the U's appearing in eq. @ are just the 
U2 's from eq. (||); that is, to leading order in the weak 
field, the plane waves in the forward light cone are just 
gauge rotated by the strong field. 

The squared amplitudes now have to be averaged over 
the density configurations. The radiation number distri- 
bution ( |l2"| ) depends on the integrated color charge den- 
sities of the sources, 



d 2 bdy 



r P (y) 



(16) 



(17) 



ill' 



Xp(y) = I dy' jUp(y') , 



Xa(v) = / dy' n 2 A (y') , 

VA 



(13) 



where y p — —yA denote the beam rapidities. Due to the 
nuclear enhancement, in the central region \A 3> Xp- Up 
to some proportionality constant which is of no relevance 
here, the integrated color charge densities Xp,a times g 
equal the saturation scales Q 2 A of the sources. 

A high-resolution probe with p t 3> Qa{v) resolves the 
gluons in the target nucleus, i.e. it appears dilute at that 
scale, while a low-resolution source with p t < Qa (y) sees 
a dense, nearly "black" target without structure (see the 
expressions for the gluon-nucleus cross section in the var- 
ious regimes in |lq]). In particular, since Qa{v) ^> Qp{y) 



Thus, fluctuations in the pt-integrated multiplicity per 
unit rapidity, dN/dy, should be dominated by fluctua- 
tions of the saturation momentum of the proton. Fluctu- 
ations of Qa do occur but their effect on dN/dy is loga- 
rithmically suppressed, and so the two quantities are ap- 
proximately uncorrelated. This can be seen also from the 
fact that for a dense nuclear target, the gluon-nucleus |l(| 
and the quark-nucleus |ll| cross sections are essentially 
given by the nuclear area ttR\ , which does not vary; once 
the target has reached the "black-body limit" , cross sec- 
tions are geometrical. Rather, high- multiplicity events 
originate from larger than average density of partons in 
the dilute projectile proton. Note that fluctuations of 
the inelastic pA cross section have been considered pre- 
viously, e.g. ref. |[{|. What is new here is that the dis- 
tribution of inelastic cross sections is related to the dis- 
tribution of correlation lengths exhibited by the effective 
action ([!]), and that the transverse momentum distribu- 
tions of produced particles changes systematically with 
multiplicity, see below. 



2 See also the results of ref. J17]], where gluon bremsstrahlung 
off nuclear targets was obtainea by resumming a Glauber mul- 
tiple collision series in the dipole model, rather than solving 
classical Yang Mills equations. 
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Thus, to first approximation an average over many 
events within one event class (which is defined by dN/dy) 
corresponds to an average^ over Qa at fixed Q p : 



dN 



d^p t dy ifi P 



dfiw(v? A ) / Vppe-Sb*-'!® 



x / Vp A e 



-S[PA,H A \ . 



■tr (\ai 



2 



(2tt) 2 

Vp A e - S lP A >^ 



(2tt) : 



■tr (|oi 



«2 



(18) 



Note the difference to fluctuations of the multiplicity at 
high transverse momentum above some (large) fixed cut- 
off pq . The integral of the distribution ( |14| ) from binary 
hard scatterings over pt > po (with po > Qa) is equally 
sensitive to fluctuations of both the projectile density as 
well as the target density |L4|. 
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FIG. 2. Schematic distribution of the pt distribution of pro- 
duced gluons in a given class of events, as determined by the 
saturation momentum of the proton Q p (or, equivalently, by 
dN/dy). Q A denotes the average saturation momentum of 
the nucleus. 

The transverse momentum distribution within one 
event class therefore has the shape depicted in Fig. ||. 
At very small transverse momentum pt < Q p , the distri- 
bution appears rather flat [jl3| . Here, both the field of the 
proton as well as that of the nucleus are in the nonlinear 
"saturation" regime. As discussed above, the width of 
this region is proportional to the square root of dN/dy. 



It is represented here as an integral over a continuous dis- 
tribution of curvatures but of course w(y?) could just be a 
sum over delta-functions if one has a potential with several 
discrete minima. 



For the highest multiplicity classes, when dN/dy is about 
200 1 / 3 ~ 6 times larger than on average, this region of 
transverse momenta could extend all the way up to Q A . 
In turn, for low multiplicity events it will shrink. 

Above Q p , up to about the average nuclear satura- 
tion momentum Q A the gluon distribution drops approx- 
imately like This is the regime where the proton 
field is weak but that of the nucleus still is "saturated" . 
It can extend down to small transverse momentum for 
the lowest multiplicity class, or could be "squeezed" com- 
pletely at very large dN/dy. 

For a closer look at the pt distributions in different 
multiplicity classes it is useful to divide the distributions 
from each class by that corresponding to the highest mul- 
tiplicity bin, forming the ratio 



R(Pf,Q 2 p ) 



d P Uy [Qp) 



dN 
dpi dy 



(19) 



This should be flat at small transverse momenta, 



Pi 



< 



Q p . Above Q p it drops since dN/dpf dy(Qp) de- 
creases while dN / dp 2 dy(Q 2 max ) is approximately con- 
stant. Above Qa finally, the ratio flattens and becomes 
constant since from eq. (^) the distributions in the var- 
ious multiplicity bins become proportional. 
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FIG. 3. pt distributions in various multiplicity bins relative 
to that in the highest-multiplicity bin. 



Fig. U shows the behavior of R(p t 
Pt distributions were taken to be 



For this figure, the 



dN 
dp\dy 



Pi 



A 2 

D2 



pUQ 2 



Q 2 pQ\ 



■A 2 ) 



pKQI 



-A 2 ) 



(pt < Qp) 

(Qp <Pt < Qa) 

(p t > Qa) , 



(20) 



(21) 



(22) 
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with the cutoff A = 0.2 GeV, the nuclear saturation scale 
Qa = 1-5 GeV and (5 Pima x = Qa- These numbers may 
be reasonable for the central rapidity region of p+Au col- 
lisions at RHIC energy. For LHC they are expected to be 
bigger by about a factor exp(A(y L HC - 2/rhic)/2) ^ 1.7, 
where j/lhc — 8.8 and ?/rhic — 5.4 are the beam rapidi- 
ties at the two colliders, and the intercept A ~ 0.3 j2(J. 

The broadening of the pt distributions with increasing 
multiplicity then leads to an increase of (p t ) with dN/dy; 
at least as long as Q p -C Qa, (pt) was shown to scale 
with the square root of dN/dy |15|l(|. This may lead to 
changes of multiparticle correlations |2lj] with increasing 
dN/dy. Moreover, at midrapidity, the inelastic cross sec- 
tion is dominated by scattering of small- a; gluons and a 
nearly flavor symmetric sea of (anti-)quarks. Therefore, 
when the typical transverse momentum of produced glu- 
ons becomes larger than typical hadronic mass scales, 
the hadronic final state should as well be nearly flavor 
symmetric. It might be impossible to achieve sufficiently 
large (pt) at either RHIC or LHC to be able to neglect, 
as an example, the masses of the kaons. Nevertheless, 
it would be very interesting to sec whether multiplicity 
ratios for e.g. K ± to tt^ 1 (or T) to 7r°) do increase at all 
with the multiplicity of the event. Experimentally, this 
is not seen for pp collisions up to Tevatron energies p2[ . 
It might be related to the fact that the correlation length 
for thermal Wilson lines at confinement is large, on the 
order of the proton radius pj| . 

Finally, it may be interesting to consider the distri- 
bution of leading hadrons (in the proton fragmentation 
region) which are produced by fragmentation of large- 
x quarks from the proton |l§|], similar to Deep Inelastic 
eA scattering 24 1. Far from the beam rapidity of the 
nucleus, its saturation momentum increases by a factor 
exp(Ay) relative to that at central rapidity, while in turn 
Q p decreases by the same factor. The quarks from the 
incident proton are then scattered to rather large trans- 
verse momenta of order Qa, and the shape of the for- 
ward p t distribution should depend less on Q p . That is, 
dN/Q 2 p dp 2 t dy should approach a "limiting shape" , which 
depends rather weakly on the event multiplicity dN/dy 
(or Qp). 

In summary, experimental analysis of the multiplicity 
distributions in high-energy pA collisions, and of corre- 
lated changes of the pt distributions could provide some 
fundamental insight into the effective action for small- 
x gluons, provided that a semi-classical distribution over 
saturation momenta exists and that it is not overwhelmed 
by higher-order corrections in a s . 



Acknowledgements: I thank Jamal Jalilian-Marian, 
Mark Strikman, and Raju Vcnugopalan for useful dis- 
cussions, the DOE for support from Grant DE-AC02- 



98CH10886, and the organizers of the workshop "Coher- 
ent Effects at RHIC and LHC: Initial Conditions and 
Hard Probes", ECT*, Trento (Italy), O ctober 14-25, 
2002, http://pluto.mpi-hd.mpg.de/trcnto, where part of 
this work was presented. 



[1] 

[2] 

[3] 

[4] 
[5] 



[6] 



[7] 
[8] 

[9] 

[10] 
[11] 



[12] 
[131 



L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rept. 
100, 1 (1983); A. H. Mueller and J. Qiu, Nucl. Phys. B 
268, 427 (1986). 

L. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 
(1994), ibid. 49, 3352 (1994); ibid. 59, 094002 (1999). 
Y. V. Kovchegov, Phys. Rev. D 54, 5463 (1996); ibid. 55, 
5445 (1997); ibid. 61, 074018 (2000). 
I. Balitsky, Nucl. Phys. B 463, 99 (1996). 
J. Jalilian-Marian, A. Kovner, L. McLerran and 
H. Weigert, Phys. Rev. D 55, 5414 (1997); J. Jalilian- 
Marian, A. Kovner, A. Leonidov and H. Weigert, Nucl. 
Phys. B 504, 415 (1997); Phys. Rev. D 59, 014014 (1999); 
ibid. 59, 034007 (1999) [Erratum-ibid. 59, 099903 (1999)]; 
J. Jalilian-Marian, A. Kovner and H. Weigert, ibid. 59, 
014015 (1999); A. Kovner and J. G. Milhano, ibid. 61, 
014012 (2000); A. Kovner, J. G. Milhano and H. Weigert, 
ibid. 62, 114005 (2000). 

E. Iancu, A. Leonidov and L. D. McLerran, Phys. Lett. B 
510, 133 (2001); Nucl. Phys. A 692, 583 (2001); E. Fer- 
reiro, E. Iancu, A. L eonidov and L. McLerr an, ibid. 703, 
489 (2002); E. Iancu, |arXiv:hep-ph/0210236 . 
A. H. Mueller, Nucl. Phys. B 558, 285 (1999). 

A. Kovner and U. A. Wiedemann, Phys. Rev. D 66, 051502 
(2002); ibid. 66, 034031 (2002). 

H. J. Pirner, Phys. Lett. B 521, 279 (2001); H. J. Pirner 

and F. Yuan, Phys. Rev. D 66, 034020 (2002). 

C. S. Lam and G. Mahlon, Phys. Rev. D 64, 016004 (2001). 

B. Svetitsky and L. G. YafTe, Nucl. Phys. B 210, 
423 (1982); R. D. Pisarski, Phys. Rev. D 62, 111501 
(2000); A. Dumitr u and R. D. Pisarsk i, Phys. Lett. 
B 525, 95 (2002); |arXiv:hep-ph/0204223| ; T. R. Miller 
and M. C. Ogilvie, Phys. Lett. B 488, 313 (2000); 
P. N. Meisinger, T. R. Miller and M. C. Ogilvie, Phys. 
Rev. D 65, 034009 (2002). 

J. Hofmann, H. Stocker, U. W. Heinz, W. Scheid and 
W. Greiner, Phys. Rev. Lett. 36, 88 (1976). 



A. Krasnitz , Y. Nara and R. Venugopalan, irXiv:hep- 
ph/0209269| ; Phys. Rev. Lett. 87, 192302 (2001) 



[14] X. N. Wang and M. Gyulassy, Phys. Rev. D 45, 844 (1992). 
[15] A. Dumitru and L. McLerran, Nucl. Phys. A 700, 492 
(2002). 

[16] A. Dumitru and J. Jalilian-Marian, Phys. Lett. B 547, 15 
(2002). 

[17] Yu. Kovchegov and A.H. Mueller, Nucl. Phys. B 529, 451 
(1998); B.Z. Kopeliovich, A.V. Tarasov and A. Schafer, 
Phys. Rev. C 59, 1609 (1999). 

[18] A. Dumitru and J. Jalilian-Marian, Phys. Rev. Lett. 89, 
022301 (2002). 

[19] L. Frankfurt and M. Strikman, Phys. Rev. Lett. 66, 2289 



5 



(1991); H. Heiselberg, G. Baym, B. Blaettel, L. L. Frank- 
furt and M. Strikman, ibid. 67, 2946 (1991). 
[20] K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 
59, 014017 (1999); J. Bartels, K. Golec-Biernat and 
H. Kowalski, ibid. 66, 014001 (2002); A. H. Mueller 
and D. N. Triantafyllopoulos, N ucl. Phys. B 640, 331 
(2002); D. N. Triantafyllopoulos, |arXiv:nep-ph/020912l| ; 
A. Freund, K. Rummu kainen, H. Weigert and A. Schaefer, 



arXiv:hep-ph/0210139 



[21] Y. V. Kovc hegov and K. L. Tuchi n, Nucl. Phys. A 708, 

413 (2002); |arXiv:micl-th/0207037 , 
[22] T. Alexopoulos el at, Phys. Rev. D 48, 984 (1993). 
[23] O. Scavenius, A. Dumitru and J. T. Lenaghan, Phys. Rev. 

C 66, 034903 (2002). 
[24] L. Frankfurt, V. Guzey, M. McDermott and M. Strikman, 

Phys. Rev. Lett. 87, 192301 (2001). 



6 



